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Abstract. We present a new method to extract the beam 
shape incorporated in the pixelized map of CMB experi- 
ments. This method is based on the interplay of the ampli- 
tudes and phases of the signal and instrumental noise. By 
adding controlled white noise onto the map, the phases are 
perturbed in such a way that the beam shape manifests 
itself through the mean-squared value of the difference be- 
tween original and perturbed phases. This method is use- 
ful in extracting preliminary antenna beam shape without 
time-consuming spherical harmonic computations. 

Key words: cosmic microwave background — cosmology: 
observations — methods: statistical 



1. Introduction 

The future space mission planck will be able to mea- 
sure, with unprecedented angular resolution and sensitiv- 
ity, the cosmic microwave background (CMB) anisotropy 
and polarization at 9 frequencies in the range 30-857 
GHz. In combination with balloon-borne experiments 
such as boomerang, maxima-1, tophat and the recently 
launched space mission MAP, these observational data will 
provide a unique base for investigation of the history and 
the large-scale structure formation of the Universe. 

The accuracy of the cosmological parameter extraction 
planned for the PLANCK mission is determined by the cor- 
responding accuracy of the systematic effects. Systematic 
errors can be one of the most important sources of er- 
rors for high multipole range of the C(l) power spectrum 
(Mandolesi et al. 2000). It is well known that extraction 
of the information about cosmological parameters such as 
baryonic density f^, cold dark matter density fl c d m , Hub- 
ble constant Hq, and so on, needs additional information 
about the statistical characteristics of the measured CMB 
anisotropy signal from the sky. The pure CMB signal is 



assumed to be a realization of a random Gaussian signal 
on the sphere with power spectrum C{1). The Gaussianity 
of the CMB signal means that all its statistical properties 
are specified by its power spectrum C(l), which depend 
on I and not on the phases. 

In the framework of the CMB observations the signal 
measured by different instruments at different frequencies, 
however, displays some peculiarities in observational as 
well as in foreground manifestations. This is why a va- 
riety of the methods of the correct information extrac- 
tion from the CMB data sets are now under discussion. 
All these methods are somewhat complementary to each 
other in the future highly sensitive CMB experiments, due 
to different sensitivity of the methods to different charac- 
teristics of the signal. 

From a theoretical point of view, the power spectrum 
of the true CMB signal is independent of Fourier rings, 
meaning that it does not depend on the azimuthal number 
m. For a flat patch of the sky it corresponds to homogene- 
ity and isotropy of the signal, without angular dependency 
of the power spectrum C(k) on 9 = tan -1 (k y / 'k x ) , where 
k = (k x , k y ). In reality, the CMB signal from the sky has 
a more complicated structure reflecting some artifacts of 
the observation and different kinds of foreground contami- 
nations, which can destroy the isotropy of the power spec- 
trum. We will focus on a few important sources causing ar- 
tificial anisotropy of the map: (i) "non-Gaussianity" (inho- 
mogeneity and anisotropy) of the foregrounds in the map; 
(ii) asymmetry of the beam shape, which is now the stan- 
dard part of investigation on systematic effects; (iii) corre- 
lations of the instrumental (pixel) noise; (iv) low multipole 
modes, e.g. / ~ 2 — 10 for the whole sky (k® ~ 1, where 
is the linear size of the path of the sky) , which are statis- 
tically peculiar. Some of the above-mentioned sources of 
the C(l) anisotropies are frequency-dependent, thus their 
contributions to the maps at different frequency channels 
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of the planck are different. For example, the foregrounds 
such as dust emission, synchrotron, free-free, as well as 
bright and faint point sources, have different frequency 
dependencies and intensities at 33 and 857 GHz range, 
which definitely can manifest as some sources of errors in 
the pixel-pixel window function and in the corresponding 
correlation function of the signals. The influence of the low 
multipole modes (point (iv)) on the possible anisotropies 
of the maps in the flat sky approximation can be detected 
directly from the corresponding C(k) amplitudes of the 
power spectrum. 

This paper is mainly devoted to illustration of the 
idea about manifestation and estimation of the asym- 
metric (elliptical or more irregular) beam shapes incor- 
porated in the pixelized data, using both analysis of 
the two-dimensional spectrum and of phases of the map 



(Naselsky et al. 2000(a)). We concentrate on beam asym- 
metry estimation using simulated CMB map, which re- 
flects directly the specific of the PLANCK scan strategy, 
map making and noise level. 

There are some important issues related to the beam 
profiles of the antenna for the PLANCK Low Frequency 
Instrument (LFI) a nd High Frequency Ins trument (HFI) 
frequency channels ( Mandolcsi et al. 200d|) . For instance, 
down to the level -10 dB at LFI, the antenna shapes have 
approximately elliptical forms and peculiarities will only 
be included at higher multipoles if the level decreases 
down to -20 dB or less, according to the design of the Focal 
Plane Unit (FPU). Thus, roughly speaking, decentraliza- 
tion of the feed horns in the FPU produces the optical 
distortions of the beam shapes from the circular Gaussian 
shapes ( Burigana et al. 199S| ) . 

Beam shape influence on the accuracy of the CMB 
anisotropy C(l) extraction from the observational data is 
related to the scanning strategy and pi xelization of th e 
maps from the time-ordered data(TOD)( Wu et al. 2001 ). 
During scanning of the CMB sky the antenna beam moves 
across the sky, meaning that antenna beam is a function 
of time. After pixelization of the TOD the position of each 
pixel in the CMB map is related directly with some points 
in the time stream for which we need to obtain the infor- 
mation of the orientation of the beam and location of the 
beam center relative to each pixel. In principle, given the 
scanning strategy( which for the planck mission is under 
discussion) and the beam shapes for each frequency chan- 
nel, we would be able to model the geometrical properties 
of the pixel beam shapes and their manifestation in the 
pixel-pixel window functions incorporated in the CMB 
power spectrum C(l). However, the computational cost 
would increase dramatically due to the complicated char- 
acter of the pixel-pixel beam matrix( Maino et al. 1999[ ). 
Moreover, the scanning strategy and the instrumental 
noise combined with the systematic effects could trans- 
form the actual beam shape during the time of observa- 
tion. We then should find some peculiarities of the com- 
plicated beam shape influence on the CMB signal. 



If the response of an antenna on the measured signal is 
linear and the CMB signal and the instrumental noise are 
Gaussian and not correlated, then the information about 
the beam anisotropy obtaining by both methods (power 
spectrum analysis and phase analysis) is the same. How- 
ever in the general case the sets of encoded information 
obtained by both methods are different. Thus using both 
methods is desirable. 

The plan of this paper is as follows. In Section 2 we 
give some definitions of CMB signals and discuss the basic 
model of the PLANCK sky map. In Section 3 we introduce a 
general power spectrum and phase analysis of CMB signal 
and the concept of beam-shape extraction. In Section 4 
we describe the main idea and its analytical approach. 
The numerical results are presented in Section 5 and the 
conclusion in Section 6. 

2. Model of the PLANCK sky map 

Let us introduce the standard model of CMB experiment 
where TOD contain the information about the signal (and 
noise ) from a large numbers of the circular scans. We sup- 
pose for simplicity that all systematic errors are removed 
after a preliminary "cleaning" of the scans. In the tempo- 
ral domain the observed signal m t is the combination of 
the CMB + foreground signal d t and random instrumental 
noise nt, 



m t = dt + n t , 
where 



dt = ^ ^ B t< i m ai m Yi m (xt), 



(1) 



(2) 



with Bt.irn being the multipole expansion of the time- 
stream beam B t (x t ). In Eq. (Q) a; m is the corresponding 
multipole coefficient of the CMB + foreground signal ex- 
pansion on the sphere and Yi m is the spherical harmon- 
ics. Following Tegmark (1996) we will assume that map- 
making algorithm is linear. The signal in each pixel s p is 
then 



d t = ^ M t,pS p , 

p=0 



(3) 



where M t . p is the corresponding pointing matrix and s p 
represents the CMB plus foregrounds signals from the sky 
convolved by the pixel beam B p j m , 



oo 



/ 



1=0 m=-l 



B p j m ai m Yi m (x p ) + N p . 



N p is the instrumental noise in each pixel: 

oo m=l 
N p = Yl N lmYlm(x p ), 

1=0 m=-l 



(4) 



(5) 
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and Ni m is the coefficient of the pixel noise expansion. x p 
is a two-dimensional vector with the components (x p ,y p ) 
to denote the location on the surface of the sphere. The 
pointing matrix M tjP depends on the scanning strategy of 
the observation. Below, as a basic model, we will use the 
model of the PLANCK mission scanning strategy discussed 
by Burigana et al. (2000) with stable orientation of the 
spin axis in the Ecliptic plane and without precession of 
the spin axis. Any precession of the spin axis will cause 
additional (regular) spin-axis modulation. 

Instead of the fast rotating beam model described by 
Wu et al. (2001), we will use a stable-orientationed beam 
model during sky crossing for the spin and optical axes, 
i.e., during the rotation of the optical axis around the 
spin axis of the satellite, the orientation of the beam is 
stable with respect to the optical axis. In addition we will 
also assume that instrumental noise is non-correlated for 
a single time-ordered scan, and between scans as well. 
Definitely this model of instrumental noise is primitive, 
and needs modifications and more detailed investigation, 
but it nevertheless reflects, as shown in the next section, 
the geometrical properties of asymmetry of the beam and 
their manifestation in the pixelized maps by a given scan- 
ning strategy. Under the assumptions mentioned above we 
will use the elliptical beam shape model of Burigana et al. 
(2000). For a small part of the flat sky approximation, we 
will use the Cartesian coordinate system with the x axis 
parallel to the scanning direction and y axis perpendicular 
to it. We denote by Xt and yt the position of the center of 
the beam at the moment t. Then the beam shape can be 
written as 



B t (x - x t ) = exp 



--(RU^D-^RU 



(6) 



with 



U 



x — Xt 

y-yt 



where R is the rotation matrix which describes the orien- 
tation of the elliptical beam, 



R 



cos a sin a 
- sin a cos a 



(8) 



with a the angle between x axis and the principal axis 
of the ellipse. The D-matrix denotes the beam dispersion 
along the ellipse principal axis, expressed as 



D 



G% 

a 2 _ 



(9) 



As mentioned by Wu et al. (2001), the signal in the time 
stream is not the pixel temperature in s p itself. The ob- 
served signal in each pixel p depends on the orientation 
of the pixel beam B p (x) and the location of its center. 
For a non-symmetric spatially dependent beam, the con- 
volution of the signal with asymmetric beam immediately 



produces asymmetry coupled with the underlying signal, 
which affects the estimation of angular power spectrum. 
For beam shape above -20 dB level, where the beam shape 
is elliptical, we will use the flat sky approximation in or- 
der to demonstrate how we can estimate the beam shape 
in the phases diagram. In such an approximation we can 
describe the signal and instrumental noise on some small 
area of the sky using FFT method instead of time consum- 
ing spherical harmonic expansion. In addition we will use 
the stable pixelized beam orientation model for a small 
area of the sky to reflect the scanning strategy of the ob- 
servation. 

3. Power and phase analysis of the CMB+noise 
signal 

Under the assumption of complete sky coverage, the spher- 
ical harmonics Yi m (9 pi 4> p ) can be expressed as a product 
of the Legendre polynomials P™ and the signal can be 
written as (Burigana et al, 1998) 



E V^T^pIHco " hV- - .Y,.„) 

2=0 



21 + 1 



■in 



E 

m=l l=Max{2,m} 

[Re(c; m 4- Ni m ) cos(m^p) 
-Im(c; m + Ni m ) sw.(m4> p )}, 



pf l (cos9 p ) x 



where 
and 



(J 



(7 ) p7\cos9 p ) = J^ fT -^ Pr(co S 9 p ) 



(10) 
(11) 
(12) 



We assume for simplicity the COBE-like cubic pixelization, 
which satisfies the following symmetry properties: if 9 G 
9 P , then —9 € 9p also, and if <ft G <fi P , then tfi + ir 6 4> p . 

Using Eq. (fLG) we can introduce the phases of the 
signal on the map, by 



tan^i m = 



Im(q r , 



■Nu 



Re(c; 



N, 



In i ) 



|Y im |sin^ se 



|AW|cos*£T e ' 



(13) 



where 



= tan _1 [Im(Q m )/Re(c; m )] is the 
phase of the signal from the sky and ^f^ se = 
tan _1 [Im(A r ; m )/Re(A^i m )] is the phase of the noise. 
Provided that CMB anisotropy is a random Gaus- 
sian field, the a/ m of eq. ^| coefficient is a random 
variable with zero mean ((a; m ) = 0) and variance 
\ a lm a l' m ') — "U' 5 mm lC ( l ) where 5[ m is the standard 
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Kronecker symbol and C(l) the power spectrum. Actually, 
the realization of the random CMB signal on the sphere 
is unique, which means that in a AT(9, </>) distribution 
on the sky we have a single realization of the phases only. 
We denote the combined signal by Si m : 

Shn=Cl m +Ni m , (14) 

and the power spectrum of this unique realization by 
|<Sz m | 2 . Our task is to extract the information about the 
beam shape either using *f?i m , \Si m \ 2 or their combination. 

The manifestation of the beam shape asymmetry in 
the |5; m | 2 can be demonstrated in the following way. Let 
us consider the following function 

A*(M = <j^g>. (15) 

Qualitatively the properties of the function A 2 (Im) are 
the following. For small I the CMB signals surpass the 
noise and the value of |Si m | 2 is determined by the CMB 
signal, so for this region (A 2 (Zm)) <C 1. For larger I, \Si m \ 2 
becomes smaller. Because of the cosmic variance (fluctua- 
tions of | Sim | 2 as a spectrum of the realization of a random 
process), some l^ml 2 can become close to zero, |S/ m | 2 < e 2 
where e is a constant and f«l. This means that A 2 (Zm) 
at these points on the (/, m)-plane has maxima 

A 2 (Zm) > I (16) 

The number density of these maxima increases at 
larger I. For I of the order It, ~ O^ 1 H> 1, where 
Of, ~ FWHM of the antenna beam, the antenna affects 
the spectrum. If the antenna is asymmetric, this influence 
is also asymmetric. Thus we have the asymmetric distri- 
bution of maxima in this region of the spectral plane. For 
very large Z's, where the noise iVj m dominates the num- 
ber density of maxima, Eq. (|I^) is determined by the 
noise, and does not depend on the beam and is therefore 
isotropic. The effects are demonstrated in Fig. [j], which is 
a result of a numerical experiment. In order to show the 
beam effect from A s (k x , k y ), we add on the symmetric 
part. There are 25 contour levels between A s (k x ,k y ) = 
and 5 x 10 -3 . The coordinates represent the flat sky ap- 
proximation of the general case at the limit I,m»l (de- 
scribed in Section 5). 

Interval of modes that are sensitive to the beam asym- 
metry starts from I ~ lb and goes to infinity in ideal con- 
ditions (no pixel noise, (A^ 2 ) = 0). For the real situation 
the limit of this interval is finite and determined by the 
pixel noise, \ai m B hn \ 2 ~ \Ni m \ 2 . 

To extract the information about the beam shape from 
Fig. |j] one needs, for example, to draw the averaged iso- 
density of the distribution of maxima of A 2 satisfied to 
Eq. (fl6|). One of the possible methods of drawing this is 
described in the next section. 




Fig. 1. Shade-filled contour map of A s (k x ,k y ) = 
(|A r k| 2 ) 1 ^ 2 /|S'k|( in logarithmic scale). The x and y axis 
are k x and k y axis, respectively. The detail of the simula- 
tion is described in Section 5. For presentation reason( to 
show the relevant part of the Fourier ring of A s (k x , k y )), 
we take 25 contour levels between and 5 x 10 -3 . 

4. Phase analysis and controlled noise as a probe 
of the antenna beam shape 

In this section we will show how to estimate the antenna 
beam shape using the information contained in the phase 
distribution of the signal in the map, using a single realiza- 
tion of the phases of all (Im) modes. After the description 
of this phase method it will be clear that, in the sim- 
ple case of linear response of the antenna, Gaussian CMB 
signal and noise, the phase method is equivalent to the 
power spectrum method described above. However, as we 
mentioned in the introduction, in the general case these 
methods give different sets of information. 

For the phase analysis we will "perturb" the phases by 
adding controlled white noise into the map. We consider 
an ensemble M ( M ^> 1), where each element of the en- 
semble consists of the same realization of the CMB signal 
and pixel noise plus a random realization of a controlled 
white noise W p with the variance a^y = const and random 
phases for each (Im) mode, 

s™ = s p + W p . (17) 

We calculate the ensemble average of the squared differ- 
ence between phases ^f m of the noise-added realization 
and the initial one, ^i m - 

A 2 (lm) = {(*Z-*im) 2 )\M (18) 

The function A 2 (Zm) from Eq. ( |l8|) is considered sepa- 
rately for different values of the variance of the controlled 
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Fig. 2. The shade-filled map of the integral( Eq. (p0|)) as a 
function of two variables: phase along the horizontal 
axis and p = \Si m \/ ^/(|W; m | 2 ) along the vertical axis. The 
intensity of the gray color corresponds to 20 gray layers 
from to 7r 2 . 

noise er 2 ^ in the range from crj^ <C cr 2 N , up to the level 
of pixel noise er 2 ^ ~ ajj. We will show that the function 
A 2 (lm) reflects all asymmetric peculiarities of the initial 
signal s p . The results of such a kind of phase analysis are 
tested numerically and are presented in the next section. 
Here we describe the analytical approach to the analysis 
of the phase mixing effect to demonstrate how it is possi- 
ble to reconstruct the antenna beam shape from the phase 
distributions in (Zm)-plane. 

For the s^f signal the definition of the combined phases 



lm 



for each (lm) mode is similar to Eq 
\S lm \sm^ lm + \Wi 



\Sim\cOsVim + \Wl T 



|sin< t 



(19) 



where Si m is the multipole expansion of the combined sig- 
nal(CMB <g> BEAM + PIXEL NOISE) at the mode (lm), 
the corresponding phase, Wj m the controlled noise 
expansion, and its phase. The analytical expression 
for the function A 2 (7m) can be written in the following 
way 



A 2 (lm) = — / if 2 dip 



-p 2 /2 



-pcosip 



1 + $ 



p COS, if 

~7T~ 



where 



| Sir. 



and 

$(a-) = A f e - t2 dt. 

V7T JO 



(20) 



(21) 



(22) 



The integral in Eq. (^) has been tabulated as a function 
of two variables */,„ and p and the result is presented in 
Fig. ||. We will use this tabulation for the numerical ex- 
periment in the next section. In this section we will obtain 
simple analytical asymptotics of Eq. (|l8|). 
From Eq. (j20|) one can find the difference 



tan* 



tan *; 



*im) 



|Wi m |cos*£ 

| Sim | COS *!„ 



cos * 



(23) 



lm 



Let us determine the mean-squared value A 2 an (Zm) of the 



difference between tan and tan 



using Eq. 



,(lm) 



1 

2^ 



One finds that 



A t 2 an (M 



tan* 



(24) 



A 2 (Zm)/cos 4 * 



if A (lm) <C 1 and after integration we obtain 



2«L 



(1 



Hm 



a D 3 COs2 ^lr, 



1 - 



-^cos2*, m 



where a im = (1 - \A - 0i m )/0lm, and 



! in 



\s h 



I 2 ) i 

— COS*; m < -. 



In the asymptotic \/3i m \ « 1 we have ai r< 
and from Eq. (p3) and Eq. fl25|) we get 



At 



A 2 (lm) 



(\W h 



2\S lr 



1 - 



3(\W lr , 



4|S im | 2 cos 2 * 



■ cos 2* 



lm 



lm 



(25) 

(26) 
< 1 

.(27) 



For our approximation we can neglect the second term in 
the brackets of Eq. (E7fl and 



A 2 (lm) 



(\Wlr. 



(28) 



Properties of (IW 7 /™,! 2 ) are the same for any values of I 
(and m). Properties of l/2|S; m | 2 were discussed in the 
previous section when we discussed A 2 (lm). It is obvious 
that qualitatively the behaviors of the function A 2 (lm) 
are the same as A 2 (/m). But there is one principal differ- 
ence between the qualitative and the full correct descrip- 
tion of the asymptotic |Sz m | 2 — > 0. Namely, if |5; m | 2 — > 
then (\Wi m \ 2 )/\Si m \ 2 > 1, and the asymptotic Eq. @ is 
not valid. Under the condition \ Wi m \ 2 |>Sz m | 2 we have 
and from Eq. ( ^8[) and Eq. (|l^) we reach 



lm 



lm 



-in 



A 2 (/m)~ — -TT^l m +^i 



(29) 



In fact the controlled white noise is some kind of aver- 
aging factor in our phase analysis. An analogous method 
may be used for the power spectrum analysis of the beam 
asymmetry. 
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5. The flat sky approximation and numerical 
results 

The general properties of |S'/ m | 2 allow us to introduce more 
convenient and faster flat sky analysis, which is specifically 
useful for the antenna beam shape estimation. This model 
reflects the scanning strategy at present specified for the 
PLANCK mission, when, for a small part of the sky far from 
the North and South poles, the model of the stable and 
fixed beam orientation is adequate, without rotation and 
multi-crossing scans. The implementation of FFT signifi- 
cantly decreases the computational cost, which, for time 
consuming spherical harmonic analysis, is a major issue. 

Using definitions of the signals and noises from the 
previous section, we define S s k y (k) to be the Fourier com- 
ponent of the signal from the sky measured by antenna 
and N(k) that of the noise, which does not depend on the 
beam properties. The Fourier modulus and phases of the 
combined signal are defined as follows, 



S(k) = S sky (k) + N(k) = y/Cpj exp[i$(k)]. 



(30) 



The function C(k) is the power spectrum of the combined 
signal (observational data and noise) and <S>(k) describes 
the phases. S s k y (k) represents the CMB signal convolved 
with the beam. If the signal from the sky does not con- 
tain non-Gaussian foreground components (or its manifes- 
tation is suppressed at the "cosmologically clean " chan- 
nels from 70 to 217 GHz flMandolesi et al. 2000D ) and if 
the noise N(k) is uncorrelated white noise, the distri- 
bution of the phases ^(k) for different k must be ran- 
dom and uniform. If there are some contaminations from 
the foregrounds and(or) pixel noise is of non-Gaussian 
nature, the phase distribution over the k range can be 
weakly correlated, and can be tested by phase diagram 
method ( |Colcs and Chiang 2000| ). Such kind of uncertain- 
ties of the statistical properties of the signal C (k) are im- 
portant for separation of the pure CMB signal and noise 
for complicated characters of the beam shape. 

We start out with a squared Gaussian random map 
with the power spectrum from the angular power spec- 
trum of the ACDM model from Lee et al. (2001). 
The map simulates a 25.6° x 25.6° square realization 
of the CMB temperature fluctuations with pixel size 3 
arcmin and periodic bou ndary conditions (PBC) (see 
Bond and Efstathiou 1987 and remarks in Section 6). We 
then add pixel noise, cr noise ~ 6 x 10 -6 , after convolv- 
ing the map with an elliptical beam with long-axis and 
short-axis FWHM 12 and 9 arcmin, respectively. After 
generation of the map, which models the HFI 100 GHz 
frequency channel, we sum M = 10 2 realizations of the 
controlled noise with variance W 2 close to the pixel noise 
variance (Eq. @ and W 2 ~ (N 2 (k))) and calculate the 
mean squared difference between phase of the signal and 
that of the signal plus controlled noise in the flat sky ap- 
proximation: 




Fig. 3. The colored phase map for the A(k) function. The 
x and y axis are k x and k y axis, respectively. For presenta- 
tion reason, we show the symmetric part( Fourier ring) in 
the Fourier domain. The size of the phase map is ranged 
from k to k ( |k| = 128), where 27r|k|~ 1 ~ 12 arcmin. 

Analogously we can calculate A s (k) for the flat ap- 
proximation. In Fig. ^] we display the function A s (k), it 
is rather difficult, however, to draw the averaged con- 
tour lines for A s (k) directly from this map. It is eas- 
ier to do this for the map of the A(fc) function. Fig. |^ 
displays the A(k) function usin g color representation of 
phases flColes and Chiang 200C ). In the color representa- 
tion, phases are mapped onto the color circle, so that the 
phase of each mode is represented by a color hue. Color 
red represents 0° and the primary hues are 120° apart, i.e., 
phases of 120° and 240° are represented by color green and 
blue. Cyan and other complementary colors magenta and 
yellow represents 180° and 300°, and 60°, respectively. For 
presentation effect, we add the symmetric part (Fourier 
ring) in the Fourier domain. It is clear from the color map 
of Fig. H that in the center for large-scale modes( small k), 
the values of A(k) are approaching zero, shown by color 
red. This is due to the fact that the phases of large-scale 
modes( the inner rings) are not "perturbed" by added 
noise. On small-scale modes( the outer rings), where the 
phases are dominated by pixel noise, the added controlled 
noise randomizes the resulting phases when the controlled 
noise level is chosen the same as that of pixel noise. 

In the intermediate regime, the beam shape manifests 
itself through the A(k) function. The fuzzy regions in the 
intermediate regime reflect the anisotropy of the beam 
shape. By taking average as 



A srn (k) 



(2fc s 



1)2 ^ 

k'es 



A(k 



(32) 



A 2 (k) = ([*„(k)-v]/(k)] 



where S = {k' : {\k x - k' x \ < k s , \k y - k' \ < k s }}, we 



(31) can find the mean shape of the beam. Figure f| shows the 
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-100 -50 50 100 



Fig. 4. The contour map of A sm (k) for fc s = 7. The con- 
tour level is in unit of degree. The size of the phase map 
is taken from -k to k( |k| = 128), where 27r|k|" 1 ~ 12 
arcmin. 




Fig. 5. The colored phase map of A function for non- 
PBC map, which is extracted from the central part of the 
previous 25.6° x 25.6° PBC simulated map. The x and 
y axis are k x and k y axis, respectively. To compare with 
Fig. |[ we only show half of the fc-range of Fig. |[ i.e., from 
-k to k( |k| = 64), where 27r|k| _1 ~ 12 arcmin. 



contour map of the A sm (k) function by k s — 7. From 
the contour map, the elliptic shape is estimated roughly 
to have the same ratio of the simulated beam shape, i.e., 
4/3. 

An important issue is related to the asymmetric 
beam extraction from a real map, covering a small 
patch of the sky. Previously we used periodic bound- 
ary condition (PBC) for modeling the CMB signal 
( Bond and Efstathiou 1987 ). In reality, for some square 
9x0 patch of the whole sky map, the PBC is artifi- 
cial and one can ask how sensitive this controlled noise 
method is to the deviation of the artificial phases of the 
signal from the true distribution? To answer this question 
we show in Fig. [s| the result of a numerical experiment for 
a map which was constructed in the following way: We 
generate, as described earlier, the PBC map with the size 
O p x O p , O p — 25.6°, and extracted from this map the 
inner part O np x O np with the size O np = 12.8°. We then 
apply the controlled noise method for the non-PBC map 
and compare our results of the beam extraction with the 
PBC case. Figure || shows the colored map of this A func- 
tion. In order to compare with the PBC case, only half 
of fc-range of interest is extracted from non-PBC case, as 
the beam size relative to the map is now twice of the PBC 
case. The peculiarity of the red crossing is induced by the 
non-PBC. The difference between these two models is less 
then 5%. This implies that we can apply our method for 
small real patch of the sky directly. 



6. Conclusion 

In this paper we propose the method of power and 
phase analysis to demonstrate the manifestation of the 
anisotropy of the antenna beam shape incorporated in the 
pixelized CMB anisotropy data, and how to estimate the 
shape of the beam. 

It is worth noting that the numerical realization of the 
phase diagram method described above is based on the 
flat sky approximation and illustrates the general proper- 
ties of the CMB signal extraction from the pixelized sky 
map at the high multipole limit I 3> 1. In such a case, 
the Fourier analysis reflects directly the general proper- 
ties of the phases of a signal on the sphere, which in the 
general (whole sky) limit must be defined using spherical 
harmonic analysis. But there are a few balloon experi- 
ments (boomerang, maxima, tophat) which cover rel- 
atively small areas of the sky in comparison with planck 
and the recently launched MAP missions. Moreover, for 
the beam shape extraction from MAP and PLANCK mis- 
sions it will be convenient to extract preliminary informa- 
tion about antenna beam shape without time consuming 
spherical harmonic computations. In connection with the 
PLANCK mission this approach looks promising due to oth- 
erwise high computation cost in the framework of the C(l) 
extraction program. 

As it is shown, the functions A^(k) and A 2 (k) reflect 
the general properties of the power spectrum C (k) mea- 
sured from the small patch of the sky. 

At the end of this discussion we would like to give the 
following remark. For estimation of the asymmetry of the 
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antenna beam shape using the controlled noise method, 
we need to know the limit of the beam cross- level (in dB), 
for which we can extract the peculiarities of the beam. 
According to general prediction of the beam shape prop- 
erties for the PLANCK mission it is realistic to assume that 
the ellipticity of the beam preserves up to the cross-level 

v ~ — 10 dB. For a small patch of the map we can use 
B(k) ~ exp[-/c2/( 2 ff2)_ fc 2/( 2cr 2)] and c ^ for the 

best fit ACDM cosmological model from the MAXIMA- 
1 data. Because of logarithmic dependence of the beam 
shape parameters on the C(k) and a 2 



N- 



2o-2 20-2 



C(k) 



N 



we can find the ratio cr^/tri with 3-5 % accuracy at the 
cross-level v ~ —6 dB. This method is most applicable to 
the patch of sky at low declination where the assumption 
of stable orientation of the antenna beam is satisfied. A 
more general case will be investigated in a separate paper. 
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